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ABSTRACT 

In the previous two papers, namely, Anusha & Nagendra (2011) and Anusha et al. (2011) we 
solved the polarized radiative transfer (RT) equation in multi-dimensional (multi-D) geometries, 
with partial frequency redistribution (PRD) as the scattering mechanism. We assumed Rayleigh 
scattering as the only source of linear polarization (Q/I, U /I) in both these papers. In this paper 
we extend these previous works to include the effect of weak oriented magnetic fields (Hanle 
effect) on line scattering. We generalize the technique of Stokes vector decomposition in terms of 
the irreducible spherical tensors Tq , developed in Anusha & Nagendra (2011), to the case of RT 
with Hanle effect. A fast iterative method of solution (based on the Stabilized Preconditioned 
Bi- Conjugate- Gradient technique), developed in Anusha et al. (2011), is now generalized to the 
case of RT in magnetized three-dimensional media. We use the efficient short-characteristics 
formal solution method for multi-D media, generalized appropriately to the present context. The 
main results of this paper are the following: (1) A comparison of emergent (I, Q/I,U /I) profiles 
formed in one-dimensional (ID) media, with the corresponding emergent, spatially averaged 
profiles formed in multi-D media, shows that in the spatially resolved structures, the assumption 
of ID may lead to large errors in linear polarization, especially in the line wings. (2) The multi- 
D RT in semi-infinite non-magnetic media causes a strong spatial variation of the emergent 
(Q/I,U/I) profiles, which is more pronounced in the line wings. (3) The presence of a weak 
magnetic field modifies the spatial variation of the emergent (Q/I, U/I) profiles in the line core, 
by producing significant changes in their magnitudes. 



Subject headings: 
atmosphere 



line: formation - radiative transfer - polarization - scattering- magnetic fields - Sun: 



1. INTRODUCTION 

Multi-dimensional (multi-D) radiative transfer 
(RT) is important to advance our understanding of 
the solar atmosphere. With the increase in the re- 
solving power of modern telescopes, and the com- 
puting power of supercomputers, multi-D polar- 
ized line RT is becoming a necessity, and practi- 
cally feasible. The multi-D effects manifest them- 
selves in the resolved structures on the Sun. The 
finite dimensional structures on the solar surface 
lead to inhomogeneity in the atmosphere, which 
is then no longer axi-symmetric. The presence of 
magnetic fields adds to the non-axisymmetry, in 



the microscopic scales through the Hanle effect. 
The purpose of this paper is to address the rel- 
ative importance of non-axisymmetry caused by 
geometry, and oriented magnetic fields. 

In the past decades extensive studies on line RT 
in multi-D media are done. A historical account on 
these developments is given in Anusha & Nagen- 
dra (2011, hereafter Paper I). In Paper I we pre- 
sented a method of Stokes vector decomposition, 
which helped to formulate an 'irreducible form' of 
the polarized line transfer equation in a 3D Carte- 
sian geometry. Such a formulation is advantageous 
because, the source vector and the mean intensity 
vector become angle independent in the reduced 
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basis. Also the scattering phase matrix becomes 
independent of the outgoing directions (SI). This 
property leads to several advantages in numerical 
work. It also provides a framework in which the 
transfer equation can be solved more conveniently, 
because the decomposition is applied to both the 
Stokes source vector, and the Stokes intensity vec- 
tor. In Anusha et al. (2011, hereafter Paper II), 
we focused our attention on devising fast numer- 
ical methods to solve polarized RT equation with 
partial frequency redistribution (PRD) in a two- 
dimensional (2D) geometry. In Paper I and Paper 
II we considered the case of non-magnetic reso- 
nance scattering polarization. Manso Sainz & Tru- 
jillo Bucno (1999) and Dittmann (1999) solved the 
polarized RT equation in the presence of a mag- 
netic field (Hanle effect), in multi-D media. Their 
calculations used the assumption of complete fre- 
quency redistribution (CRD) in line scattering. In 
this paper we solve the same problem, but for the 
more difficult and more realistic case of Hanle scat- 
tering with PRD. The physics of PRD scattering 
is treated using the frequency-domain based ap- 
proach developed by Bommier (1997a, b). The RT 
calculations in one-dimensional (ID) geometry, us- 
ing this approach are described in Nagendra et al. 
(2002). We extend their work to 2D and 3D ge- 
ometries. For simplicity we restrict to the case of 
angle averaged PRD functions. 

The present paper represents a generalization 
to the magnetic case, the decomposition tech- 
nique developed in Paper I. It also represents 
the generalizations to the 3D case, the Stabi- 
lized Pre-conditioned Bi-Conjugate Gradient (pre- 
BiCG-STAB) method developed in Paper II. An- 
other generalization is the use of 3D short charac- 
teristics formal solver in this paper, for the case of 
PRD. 

In Section 2 we describe the multi-D transfer 
equation in the Stokes vector basis. The decom- 
position technique as applied to the case of a mag- 
netic multi-D media is described in Section 3. In 
Section 4 we briefly describe the 3D short char- 
acteristics formal solution method. Section 5 is 
devoted to a brief description of the numerical 
method of solution. Results and discussions are 
presented in Section 6. Conclusions are given in 
Section 7. 



2. THE POLARIZED HANLE SCAT- 
TERING LINE TRANSFER EQUA- 
TION IN MULTI-D MEDIA 

In this paper we consider polarized RT in ID, 
2D and 3D media in Cartesian geometry (see Fig- 
ure 1). We assume that the ID medium is infi- 
nite in the X and Y directions but finite in the 
Z direction. For 2D, we assume that the medium 
is infinite in the X direction, but finite in the Y 
and Z directions. The 3D medium is assumed to 
be finite in all the X, Y and Z directions. We 
define the "top surface" for a ID medium to be 
the infinite XY plane passing through the point 
2 max . For a 2D medium, the top surface is defined 
to be the plane passing through the line (Y, Z max ), 
which is infinite in X direction. For a 3D medium, 
the top surface is the plane (A, Y, Z max ) which is 
finite in X and Y directions. For a given ray with 
direction SI, the polarized transfer equation in a 
multi-D medium with an oriented magnetic field 
is given by 

n-VJ(r,n,s) = -[ni(r)4>(x) + K c (r)} 
x[I(r,Sl,x)-S(r,Sl,x)}, (1) 

where I = (I, Q, U) T is the Stokes vector, with 
/, Q and U the Stokes parameters defined below. 
Following Chandrasckhar (1960), we consider an 
elliptically polarized beam of light, the vibrations 
of the electric vector of which describe an ellipse. 
If Ii and I r denote the components of the specific 
intensity of this beam of light along two mutually 
perpendicular directions I and r, in a plane (see 
Figure 2) transverse to the propagation direction, 
then we define 

I = Il+I r , 
Q = h — Ir, 
t7 = (/,-/ r )tan2x, (2) 

where \ is the angle between the direction I and 
the semi-major axis of the ellipse. Positive value of 
Q is defined to be a direction parallel to I and nega- 
tive Q to be in a direction parallel to r. The quan- 
tity r — (x, y, z) is the position vector of the ray in 
the Cartesian co-ordinate system. The unit vec- 
tor SI — (rj, 7, /i) = (sin 9 cos <p , sin 9 sin ip , cos 9) 
describes the direction cosines of the ray in the 
atmosphere, with respect to the atmospheric nor- 
mal (the Z-axis) , where 9 and ip are the polar and 
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azimuthal angles of the ray (see Figure 2). The 
quantity k; is the frequency averaged line opacity, 
(f> is the Voigt profile function and k c is the contin- 
uum opacity. Frequency is measured in reduced 
units, namely x = (v — i/o)/Ave>, where Ai/p is 
the Doppler width. The Stokes source vector in 
a two-level atom model with unpolarized ground 
level is 

S(r,n,x) 

_ Ki(r)(f)(x)Si(r,n,x) + K c (r)S c (r,x) 
m(r)<j>(x) + n c (r) 

(3) 

Here S c is the continuum source vector given by 
(5„(r),0,0) T with B v (r) being the Planck func- 
tion. The line source vector is written as 

/+oo 
dx' 
-oo 

r d jVR(x,x>,nn>,B) 

J 4?r <j>(x) 

(4) 

Here R is the Hanle redistribution matrix and 
B represents an oriented vector magnetic field, 
e = T i / (T R + T i) with Tj and Tr being the inelas- 
tic collision rate and the radiative de-excitation 
rate respectively. The thermalization parameter e 
is the rate of photon destruction by inelastic colli- 
sions. The damping parameter is computed using 
a = a R [l + (T E + T I )/T R ] where a R = T R /4nAu D 
and is the elastic collision rate. We denote 
the thermal source vector by G(r) = eB u {r) 
with B„{r) = (B„(r),0,0) T The solid angle el- 
ement dfl' = sin 0' d6' df', where 6 <G [0,ir] and 
if G [0,27r]. The transfer equation along the ray 
path takes the form 

dl(r,fl,x) 

ds 

= -K tot (r,x)[I(r,fl,x) - S(r,Sl,x)]. 

(5) 



The formal solution of Equation (5) is given by 
I(r,n,x) 

= J(r ,rj,x)cxp K tot (r - s"Sl,x)ds"^ 

+ f S(r - «'n, «, x)n tot (r - s'ft, x) 

J SQ 

xexp|— j Ktot(r — s"£l,x)ds"^ ds'. (6) 

I(r ,£l,x) is the boundary condition imposed at 
r = (x ,yo,z ). The ray path on which the for- 
mal solution is defined is shown in Figure 3. 

3. Decomposition of S and J for multi-D 
transfer in the presence of a magnetic 
field 

As already discussed in Paper I, a decomposi- 
tion of the Stokes source vector S and the inten- 
sity vector J in terms of the irreducible spherical 
tensors is necessary to simplify the problem. In 
Paper I, it was a generalization to the 3D non- 
magnetic case, of the decomposition technique for 
the ID transfer problems, developed by Frisch 
(2007, hereafter HF07). Here we extend our work 
of Paper I to include the magnetic fields. A sim- 
ilar technique, but in the Fourier space was pre- 
sented in Faurobert-Scholl (1991) and Nagendra 
et al. (1998), who solved the Hanle scattering RT 
problem in ID geometry. The solution of polar- 
ized Hanle scattering transfer equation using the 
angle averaged and angle dependent redistribu- 
tion matrices was presented in Nagendra ct al. 
(2002), where a perturbation method of solution 
was used. A Polarized Approximate Lambda Iter- 
ation method to solve similar problems, using the 
Fourier decomposition technique was presented in 
Fluri et al. (2003), but only for the case of angle 
averaged PRD. 

A general theory of PRD for the 2-level atom 
problem with Hanle scattering was developed by 
Bommier (1997a, b). It involves the construction 
of PRD matrices that describe radiative plus colli- 
sional frequency redistribution in scattering. It is 
rather difficult to use the exact redistribution ma- 
trix R in the polarized transfer equation. For con- 
venience of applications in line transfer theories, 
Bommier (1997b) proposed 3 levels of approxima- 
tions, to handle the R matrices. In approxima- 
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tion levels 2 and 3, the R matrices were factor- 
ized into products of redistribution functions of 
Hummer (1962), and the multi-polar components 
of the Hanle phase matrix. The collisions enter 
naturally in this formalism. It is shown that such 
a factorization of R can be achieved only in certain 
frequency domains in the 2-dimensional (x, x') fre- 
quency space. In this paper we refer to this way 
of writing the PRD Hanle R matrix, as the 'do- 
main based PRD'. The definition of the domains 
are given in Bommier (1997b) (see also Nagendra 
et al. 2002, 2003; Fluri et al. 2003). We use the do- 
main based PRD, but write the relevant equations 
in a form suitable for our present context (multi- 
D transfer). We recall that in the special case of 
non-magnetic scattering, the domain based PRD 
equations for R matrix naturally go to the Domke- 
Hubcny redistribution matrix (Domke & Hubeny 
1988). We start by writing Hanle phase matrix 
in the atmospheric reference frame in terms of the 
irreducible spherical tensors for polarimetry, intro- 
duced by Landi DcgPInnocenti & Landolfi (2004, 
hereafter LL04). In this formalism the (i, j)-th el- 
ement of the Hanle phase matrix is given by 

[p ff (n,n / ,B)] ij = 

KQ Q> 

(7) 

where = (1,2,3) and 

(8) 

where the d J MM , are the reduced rotation ma- 
trices given in LL04. The magnetic Hanle Tb 
parameter takes different values in different fre- 
quency domains (see Appendix B). Tq(i,^1) are 
the irreducible spherical tensors for polarimetry 
with K = 0,1,2, -K < Q < +K (sec Landi 
Degl'Innocenti & Landolfi 2004). In this paper, 
we consider only the linear polarization. There- 
fore, K = 0, 2 and Q e [-K, +K]. For the practi- 
cal use, we need to further expand the Pr matrix 
in each of the domains in terms of Tq ■ The re- 
quired domain based expansions of the PRD ma- 
trices in terms of Tq were already given in HF07, 



applicable there to the case of ID Hanle transfer. 
We present here the corresponding equations that 
are applicable to the multi-D transfer, which now 
become ip dependent (in the ID case, those phase 
matrix components were ip independent). We re- 
strict our attention in this paper to the particu- 
lar case of angle averaged redistribution functions 
(approximation level 3 of Bommier 1997b). 

The ij-th element of the redistribution ma- 
trix in the atmospheric reference frame (Bommier 
1997b) can be written as 

Rij(x, x',n, n', b) = WnTgii, n) 

KQ 

x{m(x,x')P^ u (j,n',B) 

+r m (x,x')P« m (j,n , ,B)}. (9) 

The weights Wk depend on the line under consid- 
eration (see LL04). Here rn(x,x') and rm(x,x') 
arc the angle-averaged versions of redistribution 
functions (see Hummer 1962). The quantities 
P£ n (j,Sl',B) and P£ m (j,Sl',B) take different 
forms in different frequency domains. They are 
described in Appendix B. 

Denoting Gq — SxoSQoG^r), where G(r) = 
eB v (r), we can write the i-th component of the 
thermal source vector as 

G i (r) = ^T Q K (i,n)G«(r). (10) 

KQ 

The line source vector can be decomposed as 
S itl (r, n,x) = J2 Tg(i, n)S$ tl (r, x), (11) 

KQ 

where 

3 

x^W K {rn(x,x')P^ u (j,n',B) 

3=0 

+rm(x,x')P£ m (j, n'.BjJ/^r.n'y). 

(12) 

Note that the components 5g Z (r,x) now depend 
only on the spatial variables (x,y,z), frequency 
x. The (9, ip) dependence is fully contained in 
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Tq{i,^l). These quantities are listed in LL04 
(chapter 5, Table 5.6, p. 211). Substituting Equa- 
tion (11) in Equation (6), the components of J can 
be written as 

h (r, n,x) = J2 T#(i, (1)1$ (r, fi, x), 

KQ 



(13) 



where 



I%(r,n,x) = I« a (r a ,n,x)e- T 



+ 



f 

Jo 



(l- Px )S^ c (r',x) dr' x {r'). 



(14) 



Here Iq — Io(ro,Q,x)5KoSQo are the intensity 
components at the lower boundary. The quan- 
tities Sq C — Sc{r 1 x)8Ko5Qo denote the contin- 
uum source vector components. We assume that 
Sc(r, x) = B v (r). The ratio of the line opacity to 
the total opacity is given by 



p x = Ki(r)<j)(x)/K to t{r,x). 



(15) 



The monochromatic optical depth scale is defined 
as 



7"x(x,y,z)=/ «tot(r - s"Q,,x)ds", 

J s 



(16) 



where t x is measured along a given ray determined 
by the direction J7. In Equation (14) T Ximax is 
the maximum monochromatic optical depth at fre- 
quency x, when measured along the ray. 

3.1. The irreducible transfer equation in 
multi-D geometry for the Hanle scat- 
tering problem 



Q 



as 



Let S« = p x S« t + (1 - p x )S« c . 1$ and S, 
well as the phase matrix elements PQ U (j,tt',B) 
and PQ IU (j,il',B) are all complex quantities. 
Following the method of transformation from com- 
plex to the real quantities given in HF07, we define 
the real irreducible Stokes vector X = (Iq, Iq, I 2 ' x , 
I 2 ' y , I 2 '*, ^2' y ) T ano - the real irreducible source 
vector S = (S° , S%, S 2 ^, S 2 ^, S 2 2 < x , S^f '. It 
can be shown that the X and S satisfy a transfer 



equation of the form 
1 



-VI ■ VX(r,fl,x) 



K to t(r,x) 
[I(r,n,x)-S(r,x)}, 



(17) 



where S(r,x) = p x Si(r, x) + (l-p x )S c (r, x) with 
Si(r,x) = eB v {r) 

i f+°° . , r dm 



+ 



dx' I ^-w{M^{B)r u (x,x') 



+M<£{B)r lll {x 1 x')}i>(n.')I(r,VL',x'), 



(18) 



and S c (r,x)^S c (r,a;),0,0,0,0,0) T W is a di- 
agonal matrix given by 

W = diag{W , 1^2,^2,1^2,^2,^2}. (19) 

The matrix W represents the phase matrix for the 
Rayleigh scattering, to be used in multi-D geome- 
tries. Its elements are listed in Appendix D. The 
matrices M^ m (B) in different domains are given 
in Appendix C. The formal solution now takes the 
form 

X(r,n,a;) =Z(r ,n,x)e- T — 



+ / e-^^S(r',x)dr' x (r'). 
Jo 



(20) 



Here X(ro,fl,x) is the boundary condition im- 
posed at ro. 

4. A 3D FORMAL SOLVER BASED ON 
THE SHORT CHARACTERISTICS 
APPROACH 

This section is devoted to a discussion of 3D 
short characteristics formal solver. Here we gen- 
eralize to the 3D case, the 2D short characteris- 
tics formal solver that we had used in Paper II. 
A short characteristic stencil MOP of a ray pass- 
ing through the point O, in a 3D cube is shown 
in Figure 4. The point O represents a grid point 
along the ray path. The point M (or P) represents 
an intersection of the ray with one of the bound- 
ary planes of a 3D cell. The plane of intersection 
is determined by the direction cosines of the ray. 
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The length As of the line segment MO (or OP) is 
given by 

As = Az//i, if the ray hits the XY plane, 
As = Ay/7, if the ray hits the XZ plane, 
As = Ax/77, if the ray hits the Y Z plane. 

(21) 

Here Ax, Ay and Az are incremental lengths 
(positive or negative) between two successive grid 
points on the X, Y and Z directions respec- 
tively. In the short characteristics method, the 
irreducible Stokes vector X at O is given by 

Z (r,n,x) = T M (r,Cl,x)exp[-AT M ] 
+ip M (r,Sl,x)S M {r,x) 
+ip (r,n,x)S (r,x) 
+V>p(r,0,x)<Sp(r,x), 

(22) 

where 5m,o.p are the irreducible source vectors 
at M, O and P. The quantity Im is the upwind 
irreducible Stokes vector for the point O. If M 
and P are non-grid points, then <Sm.p and Im are 
computed using a two-dimensional parabolic in- 
terpolation formula. While computing them, one 
has to ensure the monotonicity of all the 6 compo- 
nents of these vectors, through appropriate logical 
tests (see Auer & Paletou 1994). The coefficients 
depend on the optical depth increments in X, 
Y and Z directions. For a 2D geometry, these 
coefficients are given in Auer & Paletou (1994). 
Here we have used a generalized version of these 
coefficients, that are applicable to a 3D geometry. 

5. NUMERICAL METHOD OF SOLU- 
TION 

In this paper we generalize the pre-BiCG-STAB 
method described in Paper II to the case of a 3D 
geometry. The present work represents also an ex- 
tension of this technique to the case of polarized 
RT in the presence of an oriented magnetic field. 
The essential difference between the 2D and 3D al- 
gorithms is in terms of the lengths of the vectors. 
In a 2D geometry it is n p xn x x ny x nz whereas in 
a 3D geometry it is n p xn x x nx x ny x nz, where 
nx,Y,z are the number of grid points in the A, Y 
and Z directions, and n x refers to the number of 
frequency points. n p is the number of polarization 



components of the irreducible vectors. In the pres- 
ence of a magnetic field, n p — 6 in both 2D and 3D 
geometries. In non-magnetic problems, n p = 4, 6 
for 2D and 3D geometries respectively. 

5.1. The Preconditioner matrix 

A description of the preconditioner matrix that 
appears in the pre-BiCG-STAB method, is already 
given in Paper II. Here we give its functional form 
applicable to the problems considered in this pa- 
per. In Paper II a single preconditioner matrix was 
sufficient to handle the non-magnetic line transfer 
problem with PRD. The presence of magnetic field 
requires the use of domain based PRD matrices, 
for a better description of the PRD in line scatter- 
ing. The method requires preconditioner matrices 
to be defined, that are suitable for each of the fre- 
quency domains. We denote the preconditioner 
matrices by M M . 

A4« =I- Px 

x^iK^rn&x^ + Al^rui^x')}, 

(23) 

where 

A^S = / (24) 
and 

K% = j^WM^{B)^{n')I{r,W,x'). (25) 

Here X(r, Sl',x') is computed using a delta source 
vector as input. The expressions for the matrices 
M$ and Mjjj in different domains are given in 
Appendix C. The matrices are block diag- 

onal. Each block is a full matrix with respect to 
x and x' . The matrices are diagonal with 

respect to other variables. 

5.2. Computational details 

To calculate the integral in Equation (18) and 
the formal solution in Equation (22), we need 
to define quadratures for angles, frequencies and 
depths. 

For all the computations presented in this pa- 
per, Carlsson type B angular quadrature with an 
order n — 8 is used. All the results are presented in 
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this paper for damping parameter a — 1CP 3 . The 
number of frequency points required for a given 
problem depends on the value of a and the opti- 
cal thickness in the X, Y and Z directions (de- 
noted by Tx, Ty and Tz). A frequency band- 
width satisfying the conditions 4>(x max )Tx « 1, 
<f>(xmax)Ty « 1 and 0(x max )T z << 1 at the 
largest frequency point denoted by x max has been 
used. We have used a logarithmic frequency grid 
with a fine spacing in the line core region, and the 
near wings where the PRD effects are important. 
We use a logarithmic spacing in the A, Y and Z 
directions, with a fine griding near the boundaries. 
We find that with the modern solution methods 
used in the calculations give sufficiently accurate 
solutions for 5 spatial points per decade. 

Computing time depends on the number of 
angle, frequency and depth points considered in 
the calculations and also the machine used for 
computations. We use the Intel(R) Core(TM) i5 
CPU 760 at 2.8 GHz processor running an un- 
parallelized code. For the difficult test case of a 
semi-infinite 3D atmosphere the computing time 
is approximately an hour for one iteration. Even 
for this difficult test case the Pre-BiCG-STAB 
method needs just 18 iterations to reach a con- 
vergence criteria of 10~ 8 . 

6. RESULTS AND DISCUSSIONS 

In this section we present the results of compu- 
tations to illustrate broader aspects of the polar- 
ized transfer in ID, 2D and 3D media. We present 
simple test cases (which can be treated as bench- 
marks), to show the nature of these solutions. In 
all the calculations we assume the atmosphere to 
be isothermal. 

We organize our discussions in terms of two ef- 
fects. One is macroscopic in nature-namely the 
effect of RT on the Stokes profiles formed in 2D 
and 3D media. Another is microscopic in nature- 
namely the effect of an oriented weak magnetic 
field on line scattering (Hanle effect). We discuss 
how these two effects act together on the polarized 
line formation. 

6.1. The Stokes profiles formed due to res- 
onance scattering in 2D and 3D media 

A discussion on the behavior of Stokes profiles 
formed in ID media with PRD scattering can be 



found in Faurobert (1988) and Nagendra et al. 
(1999). In Paper II, the nature of profiles in a 
2D semi-infinite medium is compared with those 
formed in ID semi-infinite medium for CRD and 
PRD scattering (see Figures 8 and 9 of Paper II). 
Here we discuss the emergent, spatially averaged 
X and (I,Q/I, U/I) in 2D and 3D media for PRD 
scattering. 

Figures 5 and 6 show the frequency dependence 
of the components of emergent, spatially averaged 
X in 2D and 3D media respectively. The model 
parameters are, T x = T Y = T z = T = 2 x 10 9 , 
a = 1(T 3 , T E /T R = KT 4 , Tj/Tr = KT 4 , 
k-c/ki = 10 -7 , and \i = 0.11. Our choice of col- 
lisional parameters represent a situation in which 
Hi type scattering dominates. Different curves in 
each panel represent different radiation azimuths 
ipi(i = 1,12)=60, 45, 30°, 300°, 315°, 330°, 120°, 
135°, 150°, 240°, 225°, 210°. 

Iq is the largest of all the components. For the 
chosen model parameters, all the other non-zero 
components are of the same order of magnitude. 
The components I^' x and I^' y are zero in a 2D 
geometry due to symmetry reasons (see Appendix 
B of Paper II for a proof). 

The (p dependence of the X comes from the ip 
dependence of the scattering phase matrix (ty) el- 
ements. The spatial distribution of X, on the top 
surface depends sensitively on the monochromatic 
optical depths for the ray at these spatial points. 
This is a transfer effect within the medium, for 
the chosen ray direction. In the line core frequen- 
cies (x < 3), the monochromatic optical depths 
are larger, resulting in a relatively uniform spa- 
tial distribution of X on the top surface. The ip 
dependence appears as either symmetric or anti- 
symmetric with respect to the X-axis from which 
ip is measured. Thus the spatial averaging leads 
to a weak dependence of X on the azimuth an- 
gle ip. When the averaging is performed over sign 
changing quantities like the polarization compo- 
nents, it leads to cancellation, resulting in vanish- 
ing of these components. 

The <p dependence of X in the line wings can be 
understood by considering the action of the first 
column elements of the ty matrix on I® , which is 
the largest among all the components. The ele- 
ments of ty matrix are listed in Appendix D. 1$ 
is independent of p because it is controlled by the 
element tyu which takes a constant value unity. 
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Similarly Iq is controlled by \&2i which is also in- 
dependent of ip. However we see a weak ip de- 
pendence of Iq in the wings, which is due to the 
coupling of the last 4 components to Iq , which are 
of equal order of magnitude as Iq, and are sen- 
sitive to the values of ip. The p dependence of 
I±' y and T^'* elements in both 2D and 3D geome- 
tries is controlled by sin p and cos2</5 functions 
appearing in ^41 and ^51 elements respectively. 
The distribution of angle points ip in Carlson B 
quadrature is such that among the 12 p values in 
the grid, ainip takes only 6 distinct values, and 
cos2<p takes only 3 distinct values (see Table 1). 
The components I^' x and J 2 ' y are non-zero in 3D 
geometry unlike the 2D case. Their magnitudes 
are comparable to those of l\' y and 1%'* ■ The <p 
dependence of these components are controlled by 
cosy; and sin2y> functions appearing in ^/ 3 i and 
^61 elements. In the far wings, all the compo- 
nents of X go to their continuum values, as shown 
in the inset panels of Figures 5 and 6. In a ID 
geometry Iq reaches the value of B\ (parameter- 
ized as 1 here) in the far wings where the source 
function is dominated by B\. This is because 
of the fact that the formal solution with B\ as 
source function along a given ray leads to terms of 
the form B\[l — exp (— T Xtinax )]. In ID medium 
TV max = Tutot/fJ-- This implies that for semi- 
infinite ID medium, exp (— T^ max ) = so that 
Iq = B\ in the far wings. However in semi-infinite 
2D and 3D media the distances traveled by the 
rays in a given direction at different spatial points 
on the top surface are not always the same and 
therefore exp (— T x ,max) is not always zero unlike 
the ID case. Further the radiation drops sharply 
near the edges due to finiteness of the boundaries. 
Therefore when we perform spatial averaging of 
emergent Iq over such different spatial points on 
the top surface of a 2D medium (which is actually 
a line), Iq will take a value smaller than B\. For 
a similar reason (averaging over a plane) the value 
of Iq in the far wings in a 3D medium becomes 
even smaller than the value in a 2D medium. All 
other components reach zero in the far wings be- 
cause the radiation is unpolarizcd in the far wings 
(because of an unpolarized continuum). 

The way in which the components of X depend 
on ip is different in 2D and 3D geometries (compare 
Figures 5 and 6). This is a direct effect of spatial 
averaging. In a 2D medium, spatial averaging of 



the profiles is performed over the line (Y, Z msJC ) 
marked in Figure 1, whereas in a 3D medium the 
averaging is performed over the plane {X, Y, Z max ) 
marked in Figure 1. The 2D spatial averaging ac- 
tually samples only a part of the plane considered 
for averaging in a 3D medium. Also, 2D geometry 
has an implicit assumption of front-back symme- 
try of the polarized radiation field with respect 
to the infinite X axis in the non-magnetic case, 
namely 

I(r,9,p,x) = I(r,6,n-ip,x), 
I(r, 9, it + p, x) = I(r, 9, 2ir - p, x), 
Q(r, 9, p, x) = Q(r, 9, tt - p, x), 
Q(r, 9, tt + p, x) = Q(r, 9, 2ir - p, x), 
U (r, 0, p, x) = -U (r, 6,TT-p>, x), 
U(r, 9, tt + p, x) = —U (r, 9, 2tt — ip,x), 
0e[O,7r], ¥ >e[O,7r/2]. (26) 

See Appendix B of Paper II for a proof of Equa- 
tion (26). However no such assumptions are in- 
volved in 3D geometry 

Figures 7(a), (b) and (c) show I, Q/ 1,11/ 1 
profiles in non-magnetic ID, 2D and 3D media. 
Intensity / decreases monotonically from ID to 
the 3D case, because of the leaking of radiation 
through the finite boundaries in the lateral di- 
rections which is specific to RT in 2D and 3D 
geometries. In panels (b) and (c), different curves 
represent different p values. Only one curve is 
shown in panel (a), because of the axi-symmetry 
of the radiation field in the ID medium. For the 
same reason, \U/I\m = 0. The p dependence 
of |Q//|2D,3D and |£///|2D.3D directly follow from 
those of the components of X shown in Figures 5 
and 6, and their combinations (see Appendix A 
in this paper where we list the formulae used to 
construct the Stokes vector (/, Q, U) T from the 
irreducible components of X). At the line center, 
[U/I]2D,3U ~ 0. This is because U/I is zero in 
large parts of the top surface and the positive and 
negative values of U/I at x — are nearly equally 
distributed in a narrow region near the edges. A 
spatial averaging of such a distribution leads to 
cancellation giving a net value of U/I approach- 
ing zero. This is not the case in wing frequencies 
of the U/I profile (see discussions in Section 6.3 
for spatial distribution of Q/I and U/I). 
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6.2. The Stokes profiles in 2D and 3D me- 
dia in the presence of a magnetic field 

Figures 8, 9 and 10 show all the 6 components 
of X in magnetized ID, 2D and 3D media respec- 
tively. The vector magnetic field B is represented 
by {T,6 b ,Xb) = (1,90°, 68°). The corresponding 
non-magnetic components are shown as thin solid 
lines. Different line types in Figures 9 and 10 cor- 
respond to different p. The irreducible compo- 
nents in ID geometry are cylindrically symmetri- 
cal, even when there is an oriented magnetic field. 
Therefore there is only one curve in each panel in 
Figure 8. When B — the 4 components il^f 
become zero due to axi-symmetry in ID geometry 
(Figure 8). These components take non-zero val- 
ues in the line core when B^O. The magnitudes 
of Iq and Iq monotonically decrease from ID to 
3D. In the 2D case, the 2 components which were 
zero when B = 0, take non-zero values in the line 
core, when B ^ 0. Unlike ID geometry in 2D 
and 3D geometries, a non-zero B causes the last 
4 components to become sensitive to tp. The com- 
ponents ll y in 2D and J? x and I* y in 3D remain 
almost unaffected by B. This behavior is partic- 
ular to the present choice of B. For a different 
choice of B, the behavior of the 6 components 
may differ from what is shown in these figures. 
In all the geometries, the components go to their 
non-magnetic (Rayleigh scattering) values in the 
wings, because the Hanlc effect operates only in 
the line core region. 

Figures 11(a), (b) and (c) show spatially aver- 
aged /, Q/I, U /I in ID, 2D and 3D geometries re- 
spectively. Due to the finiteness of the boundaries 
in 2D and 3D media the value of spatially averaged 
/ decreases monotonically from ID to 3D. The de- 
pendence of Q/I and U/I on ip in ID medium is 
purely due to the tp dependence coming from the 
formulae used to convert X to /, Q and U (see 
Appendix A). In 2D and 3D media, the tp depen- 
dence comes from both, the tp dependence of the 
respective components of X, and also the above 
mentioned conversion formulae. The magnitudes 
of Q/I and U/I decrease in 2D and 3D geometries 
due to the spatial averaging process. The wings of 
Q/I and U/I in ID are insensitive to tp due to the 
inherent axi-symmetry. In 2D they become more 
sensitive to tp values. Again they become weakly 
sensitive to tp in 3D geometry. These differences 
in sensitivities of Q/I, U/I to the azimuth angle 



p in 2D and 3D geometries is due to the way in 
which the spatial averaging is performed in these 
geometries (see discussions above Equation (26)). 

6.2.1. Polarization diagrams in ID and 2D me- 
dia 

In Figure 12 we show polarization diagrams (see 
e.g., Stenflo 1994), which are plots of Q/I versus 
U/I for a given value of frequency x, ray direction 
(fj,,tp), and varying the field parameters two out 
of three at a time. We take F = 1, and vary 9b 
and xb values. For the 2D case we show spatially 
averaged quantities. 

For x = 0, the shapes of closed curves (loops) in 
the polarization diagrams are the same in both ID 
and 2D cases. When compared to the loops in ID, 
the sizes of the loops in 2D are smaller by about 
1% in the magnitudes of Q/I and U/I, which is 
due to spatial averaging. 

For x = 2.5, the shapes of the the loops in 2D 
are quite different from those for ID. For e.g., the 
solid curve in panel (d) is narrower than the one in 
panel (b) which correspond to 9b = 30°. On the 
other hand, the dash-triple-dotted curve in panel 
(d) is broader than the one in panel (b), which 
correspond to 9b = 120°. The orientation of a 
given loop with respect to the vertical line (Q/I = 
0) is a measure of the sensitivity of (Q/I , U/I) to 
the field orientation 9b- The size of a loop is a 
measure of the sensitivity of (Q/I, U/I) to the field 
azimuth \b- The values of \Q/I\2D and \U/I\2D 
can be larger or smaller than \Q/I\\d and \U/I\m 
for x = 2.5. The sensitivity of the line wing (x = 
2.5) polarization to (0b, Xb) is different in ID and 
2D geometries, when compared to the sensitivity 
of line center (x = 0) polarization. This is because 
at x — we sample mainly the outermost layers 
of the semi-infinite media. At x = 2.5 we actually 
sample internal inhomogeneities of the radiation 
field in (Y, Z) directions in the 2D case, and only 
those in the Z direction, in the ID case. We have 
noticed that the spatial distribution of Q/I,U/I 
at x — is relatively more homogeneous, than at 
x = 2.5 (see figures and discussions in Section 6.3 
for spatial distribution of Q/I and U/I). 
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6.3. The spatial variation of emergent 
(Q/I,U/I) in a 3D medium 

In Figure 13 we show surface plots oi Q/I and 
U 1 1 formed in a 3D media. The region chosen for 
showing the spatial distribution is the top surface 
plane (X, Y, Z max ). 

Figures 13(a), (b) demonstrate purely the ef- 
fects of multi-D geometry on the (Q/I, U/I) pro- 
files. In Figure 13(a) Q/I shows a homogeneous 
distribution at the interiors of the top surface 
(away from the boundaries) approaching a con- 
stant value (~ —3.6%). Large parts of the top sur- 
face contribute to the negative values of Q/I and 
only a narrow region near the edges contribute to 
positive values. The magnitudes of Q/I sharply 
raise near the edges. This is due to the finite 
boundaries of the 3D medium. Maximum value 
of \Q/I\ in these figures is <~ 6%. In Figure 13(b) 
U / 1 is nearly zero at the interiors of the top sur- 
face. Near the edges, the values of U/I sharply 
raise and \U/I\ takes a maximum value of ~ 20%. 

Figures 13(c), (d) demonstrate the effects 
of magnetic field on the (Q/I, U/I) profiles. 
The magnetic field vector is represented by 
B=(r,6> B ,xs)=(l,30 o ,68 o ). The nature of ho- 
mogeneity at the interior and sharp raise near the 
edges of the 3D surface, in the values of Q/I and 
U / 1 remain similar in both the magnetic and non- 
magnetic cases. An important effect of B is to 
significantly change the values of Q/I and U/I 
with respect to their non-magnetic values. \Q/I\ 
values are slightly reduced at the interior and Q/I 
now becomes —2.3%. Near the edges \Q/I\ is sig- 
nificantly enhanced and takes a maximum value of 
15%. The interior values of \U/I\ continue to be 
nearly zero. The \U/I\ is reduced at different rates 
near different edges. Now the maximum value of 
\U/I\ is 17%. We note that in ID geometry, for 
ji = 0.11, any magnetic field configuration always 
causes a decrease in \Q/I\ and a fresh generation 
of \U/I\ with respect to the non-magnetic values. 

Figures 13(c), (f) demonstrate the effects of 
PRD on the (Q/I, U/I) profiles. For this purpose 
we have chosen a wing frequency x = 5. The spa- 
tial distribution of Q /I and U / 1 is highly inhomo- 
geneous at the wing frequencies. This effect can 
be easily seen by comparing Figure 13(a) which 
exhibits large spatial homogeneity for x = 0, with 
Figure 13(e) which exhibits large spatial inhomo- 



geneity for a; = 5. For x = 0, the optical depth 
of the medium is large and therefore the radiation 
field in the line core becomes homogeneous over 
large volumes of the cube. The spatial inhomo- 
geneity of the Q/I at x = 5 is actually caused by 
the the nature of PRD function used in our compu- 
tations (which is dominated by the m function). 
Due to the frequency coherent nature of m, the 
photons scattered in the wings get decoupled from 
the line core radiation field. As the optical depth 
of the medium in the line wings is smaller than 
in the line core, the wing radiation field becomes 
more inhomogeneous and more polarized. Same 
arguments are valid for the inhomogeneous distri- 
bution of U/I on the top surface of the 3D cube. 
This can be seen by comparing Figure 13(b) with 
Figure 13(f). We recall that under the assumption 
of CRD, the values of Q/I and U/I are zero in 
the line wings (see Figure 9 of Paper II for a com- 
parison of emergent, spatially averaged Q/I, U/I 
profiles for CRD and PRD in a multi-D medium). 
The sharp increase in magnitudes of Q/I and U/I 
near the edges is larger for x = 5 when compared 
to those for x = 0. Maximum value of \Q/I\ is 
now 10% and that of \U/I\ is 40%. 

In Figure 14 we show spatial distribution of I, 
Q/I and U/I on the top surface of two different 
kinds of 3D media. Here we have chosen B = 
which is equivalent to the choice of a vertical mag- 
netic field parallel to the Z axis (because, for this 
field geometry the Hanle effect goes to its non- 
magnetic Rayleigh scattering limit). In view of the 
possible applications, we consider a cuboid with 
T x = T Y = 2 x 10 6 , T z = 20 in the left pan- 
els (a, b, c) and a cuboid with Tx = Ty = 20, 
T z = 2 x 10 6 in the right panels (d, c, f). They 
represent respectively a sheet and a rod like struc- 
ture. For the chosen optical thickness configura- 
tions, the radiative transfer effects are mainly re- 
stricted to the line core (x < 3) for the ray emerg- 
ing from the top surface. We show the results for 
x = 3 (in the left panels) and x = 1 (in the right 
panels), the frequencies for which the magnitudes 
of Q/I and U/I reach their maximum values. 

In Figures 14(a) and (d) the intensities reach 
saturation values in the interiors of the top surface 
and drop to zero at two of the visible boundaries 
(where a boundary condition of zero intensity is 
imposed for our chosen ray emerging at the top 
surface). 
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In Figures 14(b) and (c) we see that Q/I and 
U/I take values < 1% everywhere on the top sur- 
face. The magnitude of Q/I and U/I for this 
case are relatively less than those for the semi- 
infinite 3D atmospheres (compare with Figure 13). 
This can be understood using the following argu- 
ments. We are showing the results for a ray with 
(fi, ip) = (0.11, 60°) emerging from the top surface. 
The top surface for this figure refers to tz = 
where tz is the optical depth measured inwards 
in the Z direction. Using equations given in Ap- 
pendix A we can write approximate expressions 
for Q and U at the top surface as 

Q(fi = 0.11,^5 = 60°, x) w 

^L/ 2 ( M = 0.11,^ = 60°, a:), (27) 

U(fj, = 0.11, p = 60°, x) « 
^7 1 2 ' x ( M = 0.11,^ = 60 o , a; ) 

+ ^i 2 ' y (^-0.11,^ = 60°,x). (28) 

Iq is controlled by the element \&2i = 3 cos 2 9 — 1 
(see Appendix D) which appears in the scatter- 
ing integral for Sq. The factor = 3 cos 2 — 1 
represents the probability of scattering of photons 
incident from the direction 9. For 9 = 0° or 
9 = 180° (vertical incidence) ^21 is larger in mag- 
nitude compared to the cases 9 = 90° or 9 = 270° 
(lateral incidence) . For Tz — 20 the medium is ef- 
fectively optically thin (because eTz << 1) in the 
Z direction, and therefore photons easily escape in 
this direction. Thus there are smaller number of 
photons for incidence along the vertical direction 
when compared to the effectively thick case. For 
T z = 2 x 10 6 or T z = 2 x 10 9 the medium is effec- 
tively optically thick (because eTz » 1) in the Z 
direction and therefore leaking of photons in this 
direction is reduced when compared to the case of 
Tz = 20. In this way, for large values of Tz the 
probability of photons to be incident in the verti- 
cal direction is large. Therefore, as Tz increases 
the values of Iq and hence Q /I increase. 

For the chosen line of sight, Stokes U is gener- 
ated mainly by J 2 ' x and l\ ,y ■ They are controlled 
by ^31 and ^41 elements (see Appendix D) both 
of which depend on the factor sin 29. This implies 



that \&3i and ^41 are zero for both vertical and lat- 
eral incidence of photons. These elements become 
larger when the incidence is predominant in the 
direction of 9 = 45° or 8 = 135°. Using similar ar- 
guments as above we can understand the increase 
in the values of U /I with increasing values of Tz- 
The spatial distribution of Q/I and U/I is in- 
homogeneous in both left and right panels for the 
chosen core frequencies, in contrast to the homo- 
geneous distribution observed for semi-infinite 3D 
atmospheres. The extent of inhomogenity is larger 
for the left panels which correspond to smaller 
Tz value than for the right panels. The spatial 
inhomogenity could also occur due to different 
optical thicknesses along the 3 spatial directions 
leading to different number of scatterings in the 
3 directions (unlike the case of Figure 13 where 
T x = T Y = Tz). In other words, the inhomo- 
geneities in Q/I and U/I can also be caused by a 
differential leaking of radiation in the X, Y and Z 
directions. 

7. CONCLUSIONS 

This paper is dedicated to certain extensions 
of our previous works (Paper I and Paper II) on 
polarized RT in multi-D media with PRD. 

First, we present a generalization of the Stokes 
vector decomposition technique developed in Pa- 
per I, to include the magnetic fields (Hanle effect). 

Secondly, we generalize to the magnetic 3D RT, 
the efficient iterative method called the Pre-BiCG- 
STAB developed in Paper II for the non-magnetic 
2D RT. 

Thirdly, we use the more efficient 2D and 3D 
short characteristics formal solutions, with ap- 
propriate generalizations to the present context. 
With the linear formal solver used in Paper I, prac- 
tically it is difficult to compute the solutions in 
semi-infinite media. It is not the case with the 
short characteristics former solution method. In- 
deed, the solutions presented in this paper for the 
difficult case of semi- infinite media, prove this fact. 

We present several benchmark solutions com- 
puted using the code, with all the above mentioned 
generalizations. The main results of these solu- 
tions are the following. 

The emergent (Q/I, U/I) profiles in ID media 
and the emergent, spatially averaged (Q/I, U/I) 
profiles in 2D and 3D media differ significantly, 
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both in non-magnetic and magnetic cases. The 
differences are more pronounced in the wings of 
the (Q/I, U/I) profiles. The differences between 
the emergent, spatially averaged (Q/I,U/I) pro- 
files in 2D and 3D media are negligible in Q/I, 
but noticeable in U/I. 

In the non-magnetic case, at line center, the 
spatial distribution of Q/I and U/I is homoge- 
neous in the interior of the top surface, but sharply 
raise near the edges. This is purely a multi-D ge- 
ometric effect. The presence of a magnetic field 
modifies this distribution by causing a depolariza- 
tion (decrease in the magnitude) or re-polarization 
(increase in the magnitude) of Q/I and U/I. This 
is a natural consequence of the Hanle effect. In the 
line wing frequencies, magnetic and non-magnetic 
spatial distributions look the same, as Hanle ef- 
fect is confined to the line core. However, the 
spatial distribution in the line wing frequency is 
more inhomogeneous, and the sharp raise of Q/I 
and U/I near the edges is more enhanced, as com- 
pared to those at the line center. This behavior at 
line wings is mainly due to the PRD effects. These 
characteristics are not noticeable if the CRD as- 
sumption is used in line formation studies. 

We have developed efficient techniques to solve 
polarized RT in multi-D media with PRD as the 
scattering mechanism. In future, we try to apply 
these methods to understand the linear polariza- 
tion observed in the spatially resolved structures 
on the Sun. 

We would like to thank Prof. H. Frisch for 
useful suggestions which helped to improve the 
manuscript. We thank Dr. Sampoorna for use- 
ful discussions. 
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A. EXPANSION OF STOKES PARAMETERS INTO THE IRREDUCIBLE COMPO- 
NENTS 

The Stokes parameters and the irreducible Stokes vector are related through the following expressions. 
They are already given in Frisch (2007). However we present these expressions here for an easy reference. 

I(r,ft,x) = /° + cos 2 9 - 1)J 2 

— v / 3cos0sin0(/ 2 ' x cos<p — 7 2,y sini^s) 
+^(1 - cos 2 6>)(7 2 2 ' x cos 2if - I% y sin 2(p), 

(Al) 



Q(r,ft,x) = -^=(l-cos 2 tf)/ 2 

— V / 3cos0sin0(J 2 ' x costp — 7 2,y sin ip) 

— ^^(1 + cos 2 6*)(/ 2,x cos 2(p — 7 2,y sin 293), 



(A2) 



U(r, ft, x) = v / 3sin0(/ 2 ' x sin< ) 5 + J 2 ' y cos ip) 

+V3 cos 9(l2 X sin 2ip + cos 2ip). (A3) 
The irreducible components in the above equations also depend on r, ft, x and B. 

B. THE REDISTRIBUTION MATRICES IN THE IRREDUCIBLE TENSORIAL FORM 

In this paper we use the redistribution matrices defined under the approximation level III of Bommicr 
(1997b). The expressions listed below are already given in Bommier (1997b). We give them here for the 
sake of completeness. The branching ratios (see Bommier 1997b) are given by 

0"° = r , + £> + r, - <B2 » 

with £>(°) = and = cT E , where c is a constant, taken to be 0.379 (see Faurobert-Scholl 1992). 
The Hanle T B coefficient (see Bommier 1997b) takes two different forms, namely 

r B = t' k = p^r, r B = r" = a r, (B3) 

with 

— IS 

where eB/2m e is the Larmor frequency of the electron in the magnetic field (with e and m e being the charge 
and mass of the electron). Here B is the magnetic field strength. The expressions for the redistribution 
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matrices given in Bommier (1997b) involve a cut-off frequency v c (a), which is given by the solution of the 
equation 

^ = (B5) 

and a constant z = 2\/2 + 2 coming from the angle-averaging process. 
If 

zv c {a)\x'\ - (x 2 + x' 2 ) < (z - l)v 2 (a) and 
zi> c (a)|x| — (x 2 + x' 2 ) < (z — l)v 2 (a) and 

\x'\< V2v c (a) and \x\ < V2v c (a), (B6) 
then domain 1 : 



P« lu (j,n',B) = J2{/3 {K) M« Q ,(B;T' K ) 
Q' 

-aM« Q ,(B;r")}(-l) Q 'T K Q ,(3,n>), 



= Y. M QQ>M B )^y{3M)- (B7) 

Q' 

elseif 

|x'| < v c (a) or \x\ < v c (a), (B8) 
then domain 2 : 



J2M$ Q ,(B;T' K )(-lf'T« Q ,(j,n'), 



x 

Q' 

7^(2) 



^M^/ in (B)(T^)*(j,f2'). (B9) 
Q' 



else domain 3 : 



P V(j,f2',B) = [l-a/r ) ]{[r ) -a] 

^^( B ;^)(- 1 )%(j'- n ') 



Q' 



+a£(-i)«'T5 Q ,(j,fi / )}, 



= E^W' ) .ni(S)(T^)*(7,n')- (BIO) 

cndif. If 

x(x + x') < 2v 2 c (a) and x'{x + x') < 2v 2 (a), (Bll) 



14 



then domain 4 : 



p£ u (j,n',B) 

= *Y, m qq> ^ m-i) Q 'r K Q , u, si'). 



Q' 

7 K(4) 



else domain 5 : 



= E^W'.n( B )C^)*0'. n '). (B12) 
Q' 



p£ u (j,n',B) = aY,(-i) Q 'T^ Q ,(j,n'), 

Q' 

= E^SmiWC^)^","')- (B13) 
Q' 



endif. 



The symbols A / 1qq/ )IIjIII (-B), i = 1,2,3,4,5 have different expressions in different frequency domains. 
They implicitly contain the respective branching ratios and the Hanle T parameter depending upon the 
domain. 



C. THE REDISTRIBUTION MATRICES IN THE MATRIX FORM 

We introduce the diagonal matrices 

a = aE, (CI) 

with E the identity matrix, 

P = d\ & g{^\^\^\^\^\^}, (C2) 

^ = dia s{ 1 -^)' 1 -#)' 1 -#)' 

1_ pi' 1 " pi' 1 " pi}' (C3) 

The real matrices M^(B) and M^j(B) have following expressions in different domains. 
In domain 1: 



= I PM(B,T' 2 ) - &M(B,T") . (C4) 



In domain 2: 



M$(B)=\\p-a\M(B,r' K )\. (C5) 
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In domain 3: 



- a 



M{B,Y' 2 ) + a\ 



(C6) 



In domain 4: 



m£\b) = &M(B,r"). 



(C7) 



In domain 5: 



M I f(B)=d. 



(C8) 



D. The scattering phase matrix in real form in the reduced basis 



The elements of the matrix ^ are already given in Appendix A of Paper I. However we have found that 
there were some typographical errors there. We give here the elements again, correcting those typographical 
errors. 

*12 *13 *14 *15 *16 \ 
*23 *24 *25 *26 



/ *11 
*12 



22 



^13 



9. * 16 



±*23 



±*24 



1*26 



*33 *34 *35 *36 

#34 #44 #45 #46 

#35 #45 #55 #56 

#36 #46 #56 #66 



(Dl) 
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where the distinct matrix elements are: 

*n = l; *i 2 = ^|(3cos 2 0-l); 

^13 = — — sin 20 cos ip\ ^14 — sin 20 sin (p; ^i 5 = sin 2 cos 2tp; 

*i 6 = sin 2 sin 293; * 2 2 = -(9cos 4 0- 12cos 2 + 5); 

*23 = sin 20(1 - 3 cos 20) cos ^; * 2 4 = --^ sin 20(1 - 3 cos 29) sin ip; 
4y2 4y2 

/3 /3 

*25 = — 7= sin 2 0(1 + 3 cos 2 0) cos 2ip; * 2 6 = 7= sin 2 0(1 + 3 cos 2 0) sin 2^; 

2y2 2y2 

* 33 = ^sin 2 0[(l + 2 cos 2 0) - (1 - 2 cos 2 0)cos2^]; 

3 3 
^ 34 = - sin 2 0(1 - 2cos 2 0)sin2^; * 35 = — sin 20[(3 + cos 20) cos tp - (1 - cos 20) cos 3<p] ; 

* 36 = -A sin20[(3 + cos 20) sin^ - (1 - cos 20) sin3<^]; 
3 

* 44 = -sin 2 0[(l + 2 cos 2 0) + (1 - 2 cos 2 0)cos2^]; 
3 

^ 45 = — sin20[(3 + cos20) sin^+ (1 - cos 20) sin 3^]; 
3 

^46 = yrj sin 20[(3 + cos 20) cos ip + (1 — cos 20) cos 3ip] ; 

<J 55 = A [(1 + 6 cos 2 + sin 4 + cos 4 0) + (1 - 2 cos 2 + cos 4 + sin 4 0) cos 4ip] ; 
* 56 = - 2 cos 2 + cos 4 + sin 4 0) sin 4^]; 

* 66 = A [(1 + 6 cos 2 + sin 4 + cos 4 0) - (1 - 2 cos 2 + cos 4 + sin 4 0) cos Aip] . (D2) 

The elements of the matrix ^ satisfy certain symmetry properties with respect to the main diagonal. Hence 
the number of independent elements are only 21. 
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Fig. 1. — The RT in ID, 2D and 3D geometries. The 2 max , (Y, Z ma , x ), and (X, Y, Z max ) represent respectively, 
the point, the line, and the plane on which the emergent solutions are shown in this paper. The corresponding 
atmospheric reference frame is shown in Figure 2. The points A and B marked on the 2D geometry figure 
represent an example of the spatial points where the symmetry of the polarized radiation field (Equation 26) 
is valid in a 2D medium. 
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Table 1: The 12-point Carlsson type B quadrature for the azimuth angle <p. The corresponding values of 
sin ip, cos ip, sin 2^ and cos2</? are given for the purpose of discussion. 



<Pi (in degrees) 


sin ip 


COSf 


sin 2ip 


cos 2ip 


30 


0.5 


0.866 


0.866 


0.5 


45 


0.707 


0.707 


1 





60 


0.866 


0.5 


0.866 


-0.5 


120 
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-0.5 


-0.866 


-0.5 


135 


0.707 


-0.707 


-1 
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0.5 


-0.866 


-0.866 


0.5 


210 


-0.5 


-0.866 


0.866 


0.5 


225 


-0.707 


-0.707 


1 
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-0.866 


-0.5 


0.866 


-0.5 


300 


-0.866 


0.5 


-0.866 


-0.5 


315 


-0.707 


0.707 


-1 





330 


-0.5 


0.866 


-0.866 


0.5 




Fig. 2. — The atmospheric reference frame. The angle pair (9, ip) define the outgoing ray direction. The 
magnetic field is characterized by B = (T,9b,Xb), where T is the Hanle efficiency parameter and (9b, Xb) 
defines the field direction. is the scattering angle. 
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z, 




Fig. 3. — The definition of the position vector r and the projected distances r — s'fl which appear in 
Equation 6. Tq and r are the arbitrary initial and final locations that appear in formal solution integral 
(Equation 6). 




Fig. 4. — An elemental cube, showing the transfer along a section of the ray path, called a short characteristic 
(MOP). The quantities <S, K to t at M and P, and 1m at M are computed using parabolic interpolation 
formulae as M and P are non-grid points. 
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Fig. 5. — The emergent, spatially averaged irreducible Stokes vector components formed in a non-magnetic 
2D medium. Different curves represent different values of the radiation azimuth ip. The value of /i = 0.11. 
The other model parameters are given in Section 6.1. The inset panels show the far wing behavior of X. The 
x grid for these inset panels is < x < 600. 
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Fig. 6. — Same as Figure 5 but for a 3D medium. 
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(a) ID; 5 = (b) 2D; 5=0 (c) 3D; 5=0 
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Fig. 7. — The emergent, spatially averaged (I,Q/I,U/I) in non-magnetic ID, 2D and 3D media. Different 
curves represent different values of the radiation azimuth ip. The value of fj, = 0.11. The other model 
parameters are given in Section 6.1. 



24 



ID; 5=0; B*0 



ID; B=0; B*0 
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Fig. 8. — Same as Figure 5 but for a magnetic ID medium. The vector magnetic field is represented by 
(T,$b,Xb) = (1,90°, 68°). The thin solid lines show the corresponding non-magnetic components. 
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2D; 5=0; B*0 



2D; B=0; B*0 
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Fig. 9. — Same as Figure 8 but for a 2D medium. 
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Fig. 10. — Same as Figure 8 but for a 3D medium. 
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(a) ID; B^O (b) 2D; B^O (c) 3D; B^O 
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Fig. 11. — A comparison of emergent /, Q/I and U/I profiles formed in a magnetized ID media with 
the emergent, spatially averaged /, Q/I and U/I formed in a magnetized 2D and 3D media. The model 
parameters are same as in Figure 8. 
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(a) ID; x=0 (b) ID; x=2.5 
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Fig. 12. — A comparison of the polarization diagrams in ID and 2D media for two different values of 
frequency x. In 2D, the spatially averaged quantities are shown. The magnetic held parameters are given 
by r = 1, hve values of 9b in the range 30° to 150° in steps of 30°, seventeen values of xb in the range 0° 
to 360° in steps of 22°. 5. Different line types correspond to different values of 9b- Heavy square symbol 
represents xb — 0, and as we move in the counter-clockwise direction, xb takes increasingly larger values. 
The ray direction is specified by (/x, <p)=(0.11, 60°). The line types represent different 9 Bl namely (solid, 
dotted, dashed, dot-dashed, dash-triple-dotted) =(30°, 60°, 90°, 120°, 150°). 
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Fig. 13. — The spatial distribution of (Q/I,U/I) on the top surface of a 3D medium. The Q/I and U/I arc 
plotted as a function of the grid indices of tx and Ty. The ray direction is specified by (p, <^)=(0.11, 60°). 
Panels (a) and (b) demonstrate purely the multi-D effects. Panels (c) and (d) demonstrate the magnetic 
field effects. Panels (e) and (f) demonstrate the PRD effects. See Section 6.3 for details. 
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Fig. 14. — The spatial distribution of (Q/I,U/I) on the top surface of a 3D media. The Q/I and U/I 
are plotted as a function of the grid indices of tx and Ty. The ray (viewing) direction is specified by 
(fi, <p)=(0.11, 60°). Left panels represent a sheet structure and right panels represent a rod structure when 
viewed along the ±Z direction. 
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